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t^" • abstract 

o ■ 

Covariant structure of the self-force of a particle in a general curved background has been made clear in the cases of scalar 
, electromagnetic fl, and gravitational charges ]3|Q. Namely, what we need is the part of the self-field that is non- vanishing 
off and within the past light-cone of particle's location, the so-called tail. The radiation reaction force in the absence of external 
fields is entirely contained in the tail. In this paper, we develop mathematical tools for the regularization and propose a practical 
— «^ , method to calculate the self-force of a particle orbiting a Schwarzschild black hole. 

o 
cr 

?H ! I. INTRODUCTION 

, For a particle carrying a scalar, electromagnetic or gravitational charge, the field configuration of the corresponding 
• »"H . type varies in time as it moves around a black hole. To the lowest order in the charge, the particle motion follows 
a geodesic in the black hole background in the absence of external force fields. However, a part of the time-varying 
field becomes radiation near the future null infinity or future horizon and carries the energy-momentum away from 
the system, and a part of it is scattered by the background curvature and comes back to the location of the particle. 
Hence the motion of the particle is affected in the next order. The force exerted by the back-scattered self-field is 
called the local reaction force or simply the self-force. To establish a calculational strategy of this force is our ultimate 
goal. 

It is noted that we may consider the local reaction force to consist of the two parts: The part that describes the loss 
of the energy-momentum of the particle and the other part that only contributes to the shift of conserved quantities. 
In the case of geodesic motion in the Schwarzschild background, this division is unambiguous because the geodesic 
motion is completely determined by the two conserved quantities; the energy and the ^-component of the angular 
momentum (a geodesic can be assumed to be on the equatorial plane without loss of generality). In the case of the 
Kerr background, however, it seems unclear if the two parts can be identified uniquely because of the presence of the 
conserved quatity called the Carter constant which has no explicit relation to the energy-momentum of the system. 

When we attempt to calculate the reaction force on a point charge (particle), we encounter the divergence of the 
force. Hence, it is necessary to extract out the physically meaningful finite part of the force. Since the force is a 
vector by definition with respect to a background space-time, and any vector depends on the choice of coordinates in 
a covariant manner, the finite reaction force should be given covariantly. 

The covariant structure of the reaction force was investigated in the scalar case in |Q , in the electromagnetic case 
in 0], and in the gravitational case in In these investigations, the divergent part of the force was found to 

be described solely with the local geometrical quantities, whereas the finite part that contributes to the equation of 
motion was found to be given by the tail part which is due to the curvature scattering of the self-field. 

Since the tail part depends non-locally on the geometry of the background spacetime, it is almost impossible to 
calculate it directly. However, for a certain class of spacetimes such as Schwarzschild/Kerr geometries, there is a way 



to calculate the full field generated by a point charge. Considering a field point slightly off the particle trajectory, it is 
then possible to obtain the tail part by subtracting the locally given divergent part from the full field. Thus denoting 
the field by s (f> for the scalar (s = 0), electromagnetic (s = 1) or gravitational (s = 2) case, with its spacetime indices 
suppressed, the reaction force is schematically given by 

F a (r ) = lim F a [ s cj> tail }(x) , (1.1) 

F a [ s ^ il ](x) = F a [ s ^ n ](x) - F a [ s ct> dil }(x) , (x + z{r)) , (1.2) 

where z is the orbit of the particle with the proper time r, and r is the proper time at the orbital point at which 
we calculate the force. The symbol s tal1 stands for the tail field induced by the particle which is regular in the 
coincidence limit x — > z(t), s fu11 for the full field, and s <fi dn ' for the direct part as defined in Refs. Both s fu11 

and s </> dlr diverge in the coincidence limit x — > z(t). F a [...] is a tensor operator on the field, and is defined as 

[qPjVtf (s = 0), 

F a [ s (j>] = I ePj(<p r ,8 - (f>0- 1 )V (s = 1), (1.3) 

[ -mPc? ((/)p T: S - 55/3 7 0%;5 - \(t>i&;l3 + \g~t54? erf) VV S (s = 2), 

where P a ^ = 5^ + V a V^ is the projection tensor with V a being an appropriate extension of the four velocity v a (To) 
off the orbital point. 

In practice, it is a non-trivial task to perform the subtraction of the direct part, which we call the 'Subtraction 
Problem' . In this paper we propose a method to carry out this subtraction procedure covariantly. 

It should be noted, however, that solving the subtraction problem is not enough when one deals with the gravi- 
tational case. In the scalar or electromagnetic case, the reaction force is a gauge-invariant notion. In contrast, the 
reaction force in the gravitational case does depend on the gauge choice. Therefore, one has to fix the gauge appro- 
priately and evaluate the full metric perturbation and its direct part in the same gauge before calculating the force. 
We call this the 'Gauge Problem\ which seems to be a very difficult problem to solve. We do not discuss the possible 
solution of the gauge problem in this paper, but leave it for future work. Instead, we only calculate the direct part of 
the linear gravitational force under the harmonic gauge condition. 

In this paper, as a first step, we consider the case that the background is approximated by the Schwarzschild 
blackhole and use the Boyer-Lindquist coordinates, 



ds 2 = - f 1 - — J dt 2 +{l J dr 2 + r 2 (d6 2 + sin 2 9d(/) 2 ) . (1.4) 

We use the notation that x = {t, r, 9, cf>} stands for a field point, and z(t ) — z a — {t , r , 0, 4>} for an orbital point. 

The paper is organized as follows. In Section we describe basics of the issue. In Section III, we describe our 
strategy of regularization, which we call the 'mode decomposition regularization'. In Section IV , focusing on the 
scalar case, we present our method for the mode -decomposition of the direct part of the self-force. In doing so, we 
obtain a useful mathematical formula, Eq. (|4.30 ), for the harmonic decomposition. In Section 0, summing over the 



azimuthal modes m, we compare our result with Barack and Ori [ jD^lM . We find complete agreement between the 
two. In Section N% we discuss the pros and cons and conclude the paperby pointing out some future issues. Technical 
details as well asme mode decomposition regularization of the electromagnetic and gravitational cases are differed in 
Appendices. 



II. BASICS 



There are two important issues in the regularization. One is how to implement a regularization method in Eq. fll.l| ). 
Both s (t> lull (x) and s (f> d "(x) diverge in the coincidence limit x — > Zq. Such divergent quantities ar e dif ficult to treat 



in actual calculations, particularly in numerical computations. We discuss t his i ssue in Subsection II A. The other is 



how to evaluate the direct part of the field, which we discuss in Subsection II B 



A. Infinite series expansion 



We consider the case of a particle in geodesic orbit on the Schwarzschild spacetime. Even in the Newtonian limit, 
the integration of the orbit involves an elliptic function. Thus, numerical computaions will be necessary at some stage 
of deriving the self-force. However, we have divergence to be regularized which is difficult to treat numerically. The 
idea to overcome this difficulty is to replace the divergence by an infinite series, each term of which is finite. 

Suppose we have a unique decomposition method applicable to F a [ s 4> ml ](x) , F a [ s <fi dlT ](x) and F a [ s (f> tall ](x) as 
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fj s O(^E^[ s n(x), 

A 
A 
A 



(2.1) 
(2.2) 
(2.3) 



Because of the uniqueness of the decomposition, we have 

FX4> tm }(x) = i^ full ](x) ~ fX^Kx) ■ (2.4) 

At this stage, we assume that each term of the infinite series is finite, then it is possible take the coincidence limit 
x — > zq since i 7 ^[ s tall ](a;) is guaranteed to be finite. Therefore we have 

F a (r ) = J2^is^ il }(zo)- (2.5) 

A 

This approach itself does not jus tify a numerical method in the regularization calculation. Howeve r, b ecause of the 
convergence in the infinite sum ( |2.5[ ), we can expect that the sum of a finite number of terms in ( |2.5| ) gives us an 
approximated value of the self- force F a (ro). 

However, there is a very delicate problem in this approach. The exact decomposition calculation usually needs the 
global analytic structure of the field so that we can uniquely define each term in the infinite series. On the other 
hand, the regularization scheme is derived just by the local analysis of the field (l] m. Thus the direct part is defined 
only in the local neighborhood of the particle, and we have an ambiguity in the definition of the direct part. Because 
of this ambiguity, each term in the infinite series expansion is no more unique but depends on a global extension of 
the direct part we adopt. Nevertheless, the final result of the self-force should be unique. 

In this paper, we present a decomposition method based on the spherical harmonic series expansion. Although we 
have no explicit proof for the uniqueness of the resulting regularization counter terms for the self-force, the fact that 
our result completely agrees with that of Barack and Ori p5|,f6i strongly supports the validity of our method. 



B. Direct Part of the Scalar Field 



The derivation of the direct part s (f> d "(x) is one of the main issues in the regularization calculation. The direct part 
of the scalar field is obtained by integrating the direct part of the retarded Green function with the source charge. 
Here we focus on the scalar case. The electromagnetic and gravitational cases are treated in the same manner, details 
of which are given in Appendix [b]. 

The direct part of the retarded Green function G dlr is given in a covariant manner as 

G dir (x,x') = -^-9[T,(x),x'yA(x,x')6(a(x,x')) , (2.6) 

where a(x,x') is the bi-scalar of half the squared geodesic distance, A(x,x') is the generalized van Vleck-Morette 
determinant, £(x) is an arbitrary spacelike hypersurface containing x, and 0[£(x),x'] = 1 — #[x',£(x)] is equal to 1 
when x' lies in the past of £(x) and vanishes when x' lies in the future. We summarize the basic properties of the 
bi-scalars cr(x, x') and A(x, x') in Appendix |A| 



The physical meaning of the direct part is understood by the factor 0[£(x), x']<J(<t(x, x')j in Eq. (2.6). Since c(x, x') 
describes the geodesic distance between x and x' , the direct part of the Green function becomes non-zero only when 
x' lies on the past lightcone of x. Hence the direct part describes the effect of the waves propagated directly from x' 
to x without scattered by the background curvature. 

For the actual evaluation of the direct part, several methods have been proposed. In Ref. (7]||, the direct part of 
the field is calculated by picking up a limiting contribution in the full Green function from the light cone as 

-ffull 



</> dlr (x) = lim / dTG tun (x,z(T))S(T), (2.7) 



T r ot(x)- 



where G fu11 is the retarded Green function, S(t) is the scalar charge density, and T re t(x) is the retarded time defined 
by the past light cone condition of the field point x as 

0[£(x), z( Tiet )]S(a(x, z(r rct ))) = . (2.8) 
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A number of works have been made along this approach 
when we apply this method to a general orbit. 

In Ref. p| , the direct part was evaluated using the local bi-tensor expansion technique 
direct part is expanded around the particle location as 



However, the caluculation seems rather cumbersome 



Using the bi-tensor, the 



<r r (z) = q 



1 



a ;a (x, z(T ret ))v a (T Tet ) 



(2.9) 



where the letters /i, v, ■ ■ ■ are used for the indices of the field point x, a, (3, ■ ■ ■ for the indices of the orbital point z, 
and v a (r) is the orbital four velocity at z(t). The order of the local expansion is represented by powers of y which 
is linear to the coordinate difference between the field point x and the orbital point zq. Because the full force is 
quadratically divergent, we must carry out the local bi-tensor expansion of the full field up through 0(y). 

By evaluating the local coordinate values of the relevant bi-tensors, we obtain the local expansion of the full force 
in a given coordinate system (see Appendix |X|) . As described in Ref. |6) , this may be done in a systematic manner, 
and it is possible to obtain the explicit form of the divergence for a ge ner al orbit. However, the problem is how to 
decompose it into an appropriate infinite series. This is done in Section [TV|. 



III. MODE DECOMPOSITION REGULARIZATION 



We call the regularization calculation using the spherical harmonic expansion by the mode decomposition regular- 
ization. In this section, we briefly describe the regularization procedure in this approach. 

The harmonic decomposition is defined by the analytic structure of the field on the two-sphere. However both the 
direct field and the full field have a divergence on the sphere including the particle location, the mode decomposition 
is ill-defined on that sphere. Therefore, we perform the harmonic decomposition of the direct and full fields on a 
sphere which does not include, but sufficiently close to the orbit. The steps in the mode decomposition regularization 
are as follows. 

1) We evaluate both the full field and the direct field at 

x = {t,r,6,<j>}, (3.1) 

where we do not take the coincidence limit of either f or r 

2) We decompose the full force and direct force into infinte harmonic series as 

lm 
lm 



where F a [ s $ un / dlI ](x) are expanded in terms of the spherical harmonics Y^ m {9, </>) with the coefficients dependent 
on t and rln For the direct part, the harmonic expansion is done by extending the locally defined direct force 
over to the whole two-sphere in a way that correctly reproduces the divergent behavior around the orbital point 
zq up to the finite term. 

3) We subtract the direct part from the full part in each £, m mode to obtain 

F l a m [ s ^] = (Fi m [ s ^ u ] ~ F l a m U A "]) . (3.4) 

The we take the coincidence limit x — » zq. Here we note that one can exchange the order of the procedure, i.e., 
first take the coincidence limit and then subtract, provided the mode coefficients are finite in the coincidence 
limit. 



'Rigorously speaking, the angular components of the force is expanded as Fa ~ CimYA,im where YA,i m {A = 6,(j>) are the 
vector spherical harmonics. We also note that F a [s<j>im\(x) and F a m [ B <j)](x) are different since the tensorial property of the 
operator J ? a [...] depends on the spin s of the field s cj>. 
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4) Finally, by taking the sum over the modes, we obtain the self-force as 

F (7tt) = £^ m [.0(*b)- 



(3.5) 



It should be noted that because of the divergence of the full force and direct force along a timelike orbit, the 
mode coefficients of the full force and the direct force are not uniquely defined when we take the coincidence limit 
in 3). However, the tail force is regular along the orbit and it is uniquely defined. Therefore we expect the 

non-uniqueness of the direct force does not cause a problem as long as the cincidence limit is taken consistently for 
both the full force and the direct force. 



IV. DECOMPOSITION OF THE DIRECT PART 



The advant age of using (2.9) is t hat w e have a systematic method for evaluating the direct part, which we describe 
in Subsection IV A. In Subsection IV B, we describe our method for the harmonic decomposition of the direct part. 



A. Local coordinate expansion 

Though we have the covariant form of the local bi-tensor expansion of the direct part, it is not useful for the 
derivation of the infinite series expansion of it until we evaluate it in a specific coordinate system. Here we discuss 
the method to evaluate the bi-tensors in a general regular coordinate system. 

Before we consider the local expansion in a given coordinate system, we calculate the derivative of ( [2.9] ) , and derive 
the direct part of the force with the local bi-tensor expansion using the equal-time condition]^ 



= 



— cj{x,z{t)) 



We define the extension of the four- velocity off the orbit by 

V (x) := <jaa{x, Z cq )v cq , 



(4.1) 



(4.2) 



where g a& is the parallel displacement bi-vector, z cq = z(r cq (x)), and v" q = dz a /dr| T=Tcq ( a; ). Using the formulas in 
Rcf. ffll, we have 



F a [<j> dn '](x) = qg a a (x,z cq )-^- <j a. a (x,z cq ) + ^e 2 i? a ^(z cq )wf q cr ;7 (x, z cq )v 5 cq } + 0{y) , 



^2a(x,z cq ) , 

\f ~ a ap( x i z cq)Vg q V^ q 



= 1 + ^R s ^s( z cc l )Vc q v' P {x, z cq )v2 q cr ;5 (x, z cq ) + 0(y i ) 



The bi-tensors neccesary for the evaluation of the direct force (4.3) are a(x,x) and g a &{x, x), which satisfy 



a(x,x) = ^g al3 a ;a (x,x)a. l 3(x,x) = ig Q/ V a (x, x)a.p{x, x) , 

lirn (J:a(x, x) = lim a- a (x, x) = 0, 

X — *X ' X — >x ' 

g aBt . ! i{x,x)g M {x)(j. 1 {x,x) = 0, g aa .p(x, x)g Pl (x)a.^(x, x) = 0, 
\im g a a = S a a . 



(4.3) 
(4.4) 

(4.5) 



(4.6) 
(4.7) 

(4.8) 
(4.9) 



^Actually, the equal-time condition is not essential for the evaluation of the direct part. However, one can see the dependence 
on the spin of the field in a transparent manner under the equal-time condition as we describe in Appendix [b| 
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In addition, we need the generalized van Vleck-Morette determinant, 

A(x,x) = det(—g aa (x,x)a-a/3(x,x)) 



(4.10) 



We consider the local coordinate expansion of these bi-tensors around the coincidence limit x — > x, assuming that we 
have no coordinate singularity at x. 

In the coincidence limit, the effect of the curvature is small and we know the exact forms for half the geodesic 
distance bi-scalar and the paralell displacement bi-vector in the locally Cartesian coordinates as 



a(x, x) = ^q af3 (x a - x a ){ X P - x 13 ) + 0(\x - x\ 3 ) , 
g a &{x, x) = i] aa + 0(\x - x\) . 
Therefore, in a general regular coordinate system, <j(x,x) and g a& {x,x) can be expanded as 



1 



v{x,x) = -g a p(x)y al3 + 



E • 

n=3,4,— 



-A a i a 2... an (x)y a 



Jaa(x,x) = g a& {x) + ^ -] B aa\(3 l 



/3 2 -/3' 



{x)y 



n=l,2, 



where 



(x° 



c» )(a 



x a )• 



(4.11) 

(4.12) 
(4.13) 

(4.14) 



To calculate the reaction force to a monopole particle, it is enough to know the expansion coefficients of n — 3, 4 of 
Eq. (4.12) and n = 1, 2 of Eq. (Qjj).P] For a general metric, from (4.6) and (fO), we have 



-4, 



3 

i Q /3 7 — 23(q/3,7) J 



B 



a/3|7 



■ /3|q7 •> 



i 



(4.15) 
(4.16) 
(4.17) 

(4.18) 



The explicit evaluation of these coe fficie nts in the Boyer-Lindquist coordinates is given in Appendix |A|. 

The local expansion of the force (|4.3| ) on the Boyer-Lindquist coordinates is quite tedious, though systematic. We 
implement this calculation using Maple(R), an algebraic calculation program, and have 100 pages output in the end. 
However, most of the terms make a vanishing contribution to the harmonic coefficients in the coincidence limit x — * zq. 
Below, we shall focus on the terms that are non-vanishing in the coincidence limit. 



B. Harmonic decomposition of the direct part 

Without loss of generality, we may assume that the particle is located at 9q = 7r/2, </>o = at time to. Since the 
full force is calculated in the form of the Fourier-harmonic expansion and the Fourier modes are independent of the 
spherical harmonics, we may take the field point to lie on the hypersurface t = t$ in the full force. Hence we may take 
t = to before we perform the local coordinate expansion of the direct force. That is, we consider the local coordinate 
expansion of the direct force at a point {^o, r, 9, <f\ near the particle localtion {to, rg, tt/2, 0}. 

The local expansion of the direct force on the Boyer-Lindquist coordinates can be done in such a way that it consists 
of terms of the form, 



£2n 4 - 



(4.19) 



where m, ri2, tiz, n& are non- negative integers, and 



'•In the calculation of the direct force given below, only Eq. (4.12) but Eq. (4.13) turns out to be ne cessa ry. This is a result 
of our choice of the off-world line extension of the direct force, i.e., the parallel-propagation extension (4.2). 
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with a, b and </>' defined by 



/ u \ 1/2 

f:= V2r U-cos0+-(<£-0') 2 ) , (4-20) 

i? := r - r , 9 := - | , 

(4.21) 



a 1=1 + -^^—- ^—,£ 2 R 2 , (4.22) 

2r 2 r 2 + £2 ( rQ _ 2 M) 2 v ' 

b:=C, (4.23) 
'o 

u r R, (4.24) 



,2 j_ £2 



where £ := —gudt/dr, C := g^dfi/dr, and u r := g rr dr/dr^ and is the relative angle between (0 , <ft) and (vr/2, 



There are two apparently different terms in the covariant form of the direct force given by Eq. (4.2); the first term 
in the curly brackets exhibiting the quadratic divergence, and the second term proportional to the curvature tensor 
that appears to be finite in the coincidence limit. In the local coordinate expansion, the second term will give terms 
of the form R/£ or </>/£. As shown in Section [v|, the harmonic coefficients of R/£ vanish in the coincidence limit, 
while those of are finite but they give no contribution to the final result when the infinite harmonic modes are 
summed up after the coincidence limit is taken. Hence we may focus on the first term. In passing, it is worthwhile to 
note the following fact. Since the direct force shows its dependence on the spin of the field only through this second 
term (see Appendix ||), the harmonic coefficients of the direct force, which are to be subtracted from the full force, 
will be independent of the spin of the field. 



Let us focus on the first term in the curly brackets of Eq. (4.3). Since the orbit always remains on the equatorial 
plane, the force is symmetric under the transformation 9 — » it — 9, which implies there is no term proportional to 
odd powe rs of O. Hence we only need to consider the case of ni being an even number in the general form given by 



Eq. ( 4. IE ). Then the factor <d n2 may be eliminated by expressing 6 2 in terms of £, R and <j>, and we are left with 
terms of the form, 



R" 



£2n 4 + l 



(4.25) 



Explicitly, we find 



Ft 



dir 




1 (r - 2M)£u r 1 2{r Q - 2M)££ 2 u r 
3 (r - 2M)£C i u r 



2 

F ^ = g ( 1 g ^ g Cu i 

r q [r Q (r -2M)e (r„ - 2M) 2 ? £ 3 
1 2r 2 + C 2 1 1 £ 2 1 
~2 rjj I + 2 r - 2M | 
1 (3r 2 + 4£ 2 )C 2 <\> 2 2£ 2 L 2 cf> 2 



(4.26) 



2 r-3 £ 3 r -2M£ 3 

3 (r 2 + £ 2 )£ 4 4 3 £ 2 £ 4 4 
2 r 3 , £ 5 + 2 r - 2M £ 5 ' 



(4.27) 



= This is also a result of the specific off-worldline extension chosen for the four- velocity. It is valid for the parallel-propagation 
extension, but does not hold, in general, for other extensions. 
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(4.28) 

r R i 2 , r 2\ $ 
-£u r -^ -(r +C ) — 

1 (r - 2M)Cu r 1 1 (r - 2M)(rg + AC 2 )Cu r <p 2 

+ 2 ^ J3 ^3 

3 (r - 2M)(r§ + £ 2 )£V </> 4 \ 



2 ^o 2 C 



(4.29) 



where F^ 11 = F a [cf) d "}. The absence of Fg" is because of the symmetry of the background; the orbit remains on the 
equatorial plane even under the action of the self- force. In the above, we have discarded the terms of the form R/£ 
or <j>/£. As mentioned before, and as shown in Section [v], such terms give no contribution to the final force. 

What we have to do now is to perform the harmonic decomposition of the components of the direct force given 
above. To do so, we note the following important fact. Apart from the trivial multiplicative factor of R ni which 
is independent of the spherical coordinates, the terms to be expanded in the spherical harmonics are of the form 
</>™ 3 /£ 2 ™ 4+1 , or (4> — </>')™f/£ 2 ™ 4+1 - To the order of accuracy we need (in fact only the leading order accuracy is 
sufficient; see Appendix the factor {(f> ~ 0')™ 3 ma y be eliminated by replacing it to an equivalent (^-derivative 
operator of degree acting on ^ 2n 3-2ii 4 -i j wn ich is further converted to a polynomial in m after the harmonic 
expansion of ^ 2 ™3-2n 4 -i^ Thus the only basic formula we need is the harmonic expansion of £ 2p_1 where p is an 
integer. Detailed derivation of it is given in Appendix [5| Note that, apart from the term b(<j) — <fi') 2 /2 in £ with 
respect to which we expand £ in a convergent infinite series, £ 2 p -1 is defined over the whole sphere to allow the 
straightforward harmonic decomposition. The result to the leading order in the coincidence limit a — > 1 4- is 



2vr£ £ D*- 1/2 (a)Yt m (e,<t>)Y; m (e', </>'), (4.30) 



/ c \ 2 p-i / h \p-V 2 00 1 



\V2rJ \ 2 



TP-i/a, I VTTS-P- 1/2 



e=o m=-l 

(a - l) p+1/2 , for p + | < 0, 



D tm (a)^ \ (-iy 2 P+U2 ^ r(p+l/2)r(p + n+l/2) 1 f-m 2 b\ n , ( 4 ' 31 ) 



VTTb ^ T(p + n - 1 + l/2)r(p + n + £ + 3/2) n\ \l + b 



We note that, although what we need here is only the case of an integer p, the above formula is valid for any p (except 
for the case p — —1/2). 

After the decomposition, we can take the radial coincidence limit r — > tq (followed by the angular coincidence limit 
if desired). The basic properties of the resulting mode coefficients in the coincidence limit are discussed in Section 0. 
Here we briefly explain the reason why the terms proportional to and <f>/£ give no contribution to the final result. 

The term R/£ corresponds to R times the case of p = 0, for which D^ 1 ^ 2 is finite in the limit a — ► 1 (i.e., R — ► 0). 
Hence all the coefficients vanish in the radial coincidence limit. As for </>/£, it can be replaced by (<f> — <f>')/£ which 
is equivalent to d^t; in the coincidence limit. This corresponds to the case of p = 1 multiplied by m. Hence all the 
harmonic coefficients become odd functions of m, and their sum over m for each I vanishes in the angular coincidence 
limit. As a result, the non-vanishing contribution comes only from the terms R/£, 3 , 4>/£, 3 , l/£, <f> 2 /£, 3 and <^> 4 /£ 5 . 



V. REGULARIZATION COUNTER TERMS 



In this section, we present the mode decomposition of the direct force given b y E qs. ( 4.26 ) ~ ( 4.29 ), and compare 
the resulting regularization counter terms with those obtained by Barack and Ori |16[| in their mode-sum regularization 
scheme (MSRS) §. 

Barack and Ori define the regularization counter terms as 



lim F*l r = A a L + B a + C a /L + 0{L.- 2 ) 

x=>z Q 



E 

1=0 



lim F*f -A a L-B a -C a /L 

X^Zg 



(5.1) 
(5.2) 



where F^" is the multipole ^-mode of F^", L = £ + 1/2, and A a , B a and C a are independent of L. The A a term is 
to subract the quadratic divergence, the B a term the linear divergence, and the C a term the logarithmic divergence. 
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The D a term is the remaining finite contribution of the direct force to be subtracted. As shown in Appendix |^, we 
find C a — D a — in agreement with Barack and Ori M. We also find the complete agreement of A a and B a terms 
with their results for a general geodesic orbit fnjjlfj as given below. 



The direct part of the force to be considered has the form given by Eq. ( |4.25| ) , which may be re- written as 

R ni O - ft) 113 

£2n 4 + l ' 



(5.3) 



where (774,773,774) are non-negative integers. Because the highest order of divergence is quadratic, it is sufficient to 
consider the cases n\ + 71 3 — 2n^ = —1,0 and l.|**| 



We first note that 



1 1 d ( 1 



£2« 4 +l 2n 4 - 1 r 2 + C 2 dcj) \£ 2n *- 
By using this equation recursively, we obtain 



0(2/ 4+ ) ■ (5.4) 



/\n 3 an 3 

J ^ £2n 3 -2n 4 -l _|_ ^^ri 3 -2n 4 + l^ 



In the context of the harmonic decomposition, we may replace the derivative d/dcf) by i m. Hence instead of Eq. ( |5.3| ). 
we may consider the terms of the form, 

In the above eq uatio n, we have indicated by 0(y ni+n3 ~ 2rii+1 ) the presence of correction terms of 0(y 2 ) relative to 
the original form ( |5.3| ). In terms of the regularization parameters A a , B ai C a and D a , this implies the terms with 
n\ + 773 — 2n4 = — 1 would contribute to A a and C a , and the terms with n\ + 77,3 — 2?i4 = to B a and D a , while the 
terms with rt\ + 713 — 2774 = 1 to D a . However, as shown in Appendix |^, by a general argument , we can show that 



both C a and D a vanish. Therefore we do not have to worry about the 0(y 2 ) corrections in Eq. (5.5) but may focus 
on the leading behavior of it in the coincidence limit. 



Keeping this fact in mind, we now analyze which cases of the form (|5.5[) contribute to the regularization p aram eters 



A a and B a . For this purpose, we set 771 + 773 — 2774 = q where q = —1, or 1. Then comparing Eq. (5.5) with 



Eq. (4.30), we find it is convenient to separately consider the two cases: 



(1) The case 2p = 2713 — 2774 = 773 — ni + q < —2 



In this case, the harmonic coefficients of Eq. (5.5) behave as 



~ ^ n i+2 n 3— 2n 4 + l _ j^n 3 +q+l 

Since 773 > 0, the harmonic coefficients are non- vanishing in the limit R — > only if 773 = and q = — 1. This 
means rii = 2n^ — 1 (> 0). Therefore only the terms of the form /j>2n 4 - 1^2^+1 ( n4 > give finite coefficients, 
and they contribute to A a . 

(2) The case 2p = 2713 — 2774 = 773 — ni + q > 0. 

In this case, since the harmonic coefficients of £ 2n 3-2n 4 -i are g n j^ e j n ^ ne n m n f> q, we must have n\ = 0, 
hence 773 = 2774 + q. Therefore, since D^ 1 ^ 2 is an even function of 777, the harmonic coefficents will be odd 
functions of to if q is odd, i.e., if q = —1 or 1. When the sum over to is taken, the result vanishes in the angular 
coincidence limit if q is odd because of the symmetry property of \Yi m (9,<p)\ 2 under m — » — m. Thus only the 
case of q = or 773 = 2774 remains. The corresponding terms are of the form (cf> — 0') 2 ™4^2ri 4 +i ; anc j they 
contribute to B a . 

From the above results, and noting that 4>' cx R, we obtain the equality, 

(0-0') 2n {4> 2 - 2W + 4>' 2 ) n 4> 2 



£2n+l £2n+l £2n+l 



(5.6) 



which holds in the sense of its contributions to the regularization parameters. Thus, to summarize, the non-vanishing 
contributions are from the terms either of the form /j 2 ™+ 1 /^ 2 ™+ 3 or of the form <fi 2n / £ 2n+1 , where 77 is a non-negative 
integer, and the former contributes to A a while the latter to B a . 



** A lthough we may further restrict 774 to the range < 774 < 2 from the explict form of the direct force in Eqs. (4.2£) 



(4.29), we choose not to do so because it turns out to be unnecessary in the following discussion. 



9 



A. The ,4-term 



The A-term de scrib es the quadratic divergent terms of the direct force. Thus we consider the most divergent terms 
in Eqs. flOej) - ( gjgg ), 



R 4> 

£3 and £3 • 



(5.7) 



As discussed above, the term (A/£ 3 may be replaced as {4> — </>')/£ + <f>'/£ = </>'/£■ Hence we may focus on the form 
i?/£ 3 . The essential fact is that this is odd in R. This leads to the harmonic coefficients proportional to sign(i?). 
Using the fomula (D4), we obtain 



A t = sign(i?) — 



r l-£ 2 Ao' 



-sign(-R) 



q to 



z 



r -2Ml- C 2 /rl ' 



At = 0. 



(5.8) 

(5.9) 
(5.10) 

These A-terms vanish when averaged over both limits R — > ±0. There could be correction terms of O(y ) which could 
contribute to the C and D-terms. However, as shown in Appendix they are known to be absent. 

B. The B-term 

The linearly divergent terms are described by the J3-term, which are of the form, 

(5.11) 



2n+l 



in Eqs. ( 4. 26] ) ~ ( 4. 29] ). The Legendre coefficients are given by the formula (D15). We find 



Bt = - 



(ro- 



2r Q 

(m - 2M )u 2 r 



2M)£u r ^ {2n)\ (2n + 1) 2 I> + 1/2) C 2n 

2r n ^„2 2 "(n!) 2( ' ^T(n + l) r l n+1 



n=0 



V^r(n + 1) 



(5.12) 



2r 

O rn 02" 



(2n)! ( _ 1)n (2n+l) 2 r(n+l/2) £ 



n=0 

(2n)! 



2 2 ™(n!) 2 



V / 7rr(n + 1) r. 



2n+l 

o 



2r ^ 2 2n (n!) 2 

n— 



(-(2n-l))r(n+l/2) C 2 



2n+l ' 
r 



(5.13) 



Brh — 



(r - 2M)£u r 
2T 2 



(2n)! ( _ ir (2n+l) 2 r(n + l/2) £ 



2 2 ™(n!) 2 



V5Fr(n + i) 



2n+l 
'0 



2 (2(n + l))! (2n + l) 2 I> + l/2) £ 2 («+D 

•2 / j o2(ii.+ ll//„ i 1M\2\ / 



22(n+D(( n + l)!)2' 



ViT(n + 1) 



2(ra+l) + l 



The above may be expressed in terms of the hypergeometric functions as 

_ _ {r -2M)£ur„{3 3 £ 2 
2^ F U'2 ;1; -^ 



_ (r - 2M)ug „ /3 3 ^ 
2r3 r 2 



2r 2 
z 'o 



F|II-1-^ 



Bub — 



(r - 2M)£u r 
16r| 



2 2 , 
'.2'2' A r 2 



/3 3 & 

2rf{rV Z > r 2 
9C 2 (5 5 

' 2'2' d ' r 2 



(5.14) 

(5.15) 
(5.16) 
(5.17) 



The above results for the A and £>-terms perfectly agree with the results obtained by Barack and Ori in a quite 
different fashion @,|l§. Q 



Here we give the values of B a expressed in terms of generalized hypergeometric functions, while in |l5[ these are given in 
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VI. CONCLUSION AND DISCUSSION 



Our final goal is to establish a method of calculation of the local gravitational reaction force to a point particle 
orbiting a Kerr black hole. We have pointed out in Section | that there are two problems; the 'subtraction problem' 
and the 'gauge problem'. 

In this paper, we have only discussed a possible approach to the subtraction problem. We have introduced a 
regularization method which utilizes the spherical-harmonic decomposition, and have derived the direct part of the 
self- force, which turns out to be independent of the spin s of the field under consideration. The harmonic decomposition 
of this direct part has been carried out, and the regularization counter terms for the self- force have been derived for 
a general geodesic orbit. We have found our result agrees completely with the result obtained by Barack and Ori |l(t | 
in their mode-sum regularization scheme (MSRS) |9[. 

To compare with the MSRS, we have derived the regularization counter terms which are obtained by summating 
the harmonic coefficients over m. However, when we extend our method to the Kerr background, we ma y h ave to 
carry out the regularization before taking the m-summation. In this sense, the formulas derived in Appendix Tlx where 
no summation over m is assumed, may be still useful in the Kerr case. 

It is worthwhile to point out that the gauge problem in the gravitational case seems far more serious than the 
subtraction problem. What we know at the moment is that the gravitational self-force is described by the tail part 
of the metric perturbation induced by a particle (^^]. However this is justified only in the harmonic gauge, while 
the full metric perturbation can be obtained only in the Regge- Wheeler gauge or in the radiation gauge where the 
identification of the tail part is highly non-trivial. A prescription to identify the tail part of the metric perturbation 
has been proposed in p0| , but it needs to be verified. The gauge problem for the non-radiative monopole and 
dipolc components of the metric perturbation which are not obtainable in the Teukolsky formalism seems to stand as 
additional serious obstacle. Possible resolutions for the gauge problem are under investigation. F 
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APPENDIX A: BI-TENSORS AND LOCAL COORDINATE EXPANSION 



Bi-tensors are tensors which depend on two distinct field points, say, x a and x a . For our purpose, we consider 
half the squared geodesic interval bi-scalar a(x,x), and the geodesic paralell displacement bi-vector g a a(x,x) 7 which 
satisfy 

a(x,x) = ^g a,3 a. a (x,x)(T. fj (x,x) = ^g a ^a. & (x, x)a.p(x, x) , (Al) 
\\m a - a {x , x) — lira a-„ (x, x) — 0, (A2) 

g a a-f3(x,x)g Pl (x)cr. n (x,x) = 0, g aa .p(x,x)g 01 (x)<r^(x,x) = 0, (A3) 
lim.g Q a = <S Q a . (A4) 

X — >x 

We also need the generalized van Vleck-Morcttc determinant bi-scalar, 

A(x, x) = det(-g a& {x, x)a. a( 3(x, x)) . (A5) 
We consider the local expansion of these bi-tensors around the coincidence limit x — > x. 



terms of the two complete elliptic integrals K and E. The two expressions are related by changing of variables and using the 
formulas in Jl8{ . 
"In this respect, recently some progress has been made in Jl7|. 
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calculating the local expansion of the field and its derivatives in a covariant way, the following formula are useful 

(A6) 



<T; a p(x, z) = g a p{z) - ^R a 1 f 3 S (z)a. 1 (x, z)cr. s {x, z) + 0(e 3 ) 



(T-npix, z) = -g^ a {x, z) yg a p{z) + ^R ai ps{z)a'^ {x, z)a' 5 (x 1 z) 

+ 0(6 3 ), 

gr#{x, z) = ~g^(x, z)R a 70S (z)<T' s (x, z) + 0(e 2 ) , 
^%(x, *) = -^ix, z)g^(x, z)R a M5 {z)a> 5 {x, z) + 0{e 2 ) . 
As for the calculation of A(x, x), we use the result of the covariant expansion given in Ref. @. We have 

cr ;g/3 (x, x) = -g/(x, x) (g & p{x) + ^R s ^ps{x)<r' 7 (x, x)<r ;5 (x, x) + 0{\x - x\ 3 ) 



Then, for a vacuum background, we have 



A(x,x) = 1 + 0(\x - x\ 3 ) . 



(A7) 
(A8) 

(A9) 
(A10) 



Th e local expansion of these bi-tensors by the background coordinates is derived from the formulas (4.12) and 
(4.13), a nd th e expan si on co efficiet s exp ressed in terms of the ordinary derivatives of the background metric are given 
in Eqs. (4.15), (4.16), (4.17) and ( 4.1 8| ) . In the Schwarzschild background, the non-vanishing components of these 
coefficients are 



Attr - 


- A tr t — A rt t — - 


M 




A 


3M 
(r - 2M) 2 ' 






A r gg 


= Ag r g = Agg r = 


r, 




A T( p(p 


— A^Tffj} — A(p(p r — 


= r sin 2 9 




Ag^ 


= A^g = A^g^ = 


-- r 2 sin 9 


COS0 , 


Atttt 


M 2 (r-2M) 






Attrr 


— Atrtr — A r n r - 


= At rr t — 


: A r trt = 


Atteo 


= Atete = Ague = 


- A t ggt - 


= Ag t gt = 


Att4>4> 


= A t $t<t> — A$ tt <t> 


= At^t 


= A^t - 


A 


M(Sr - M) 
r(r - 2M) 3 ' 






A rr gg 


= A r g r g = Ag rr g 


= A r gg r 


= Ag r g r 


A rr( p^ 


— A r( jj r( j) — A(p rr 


p Arcfxj 


ir A(p r( j) 



Arrtt — 
Aggtt = 



.4 



<jxj>tt 



M(4r - 5M) 
3r 3 (r-2M) 
M(r - 2M) 
3^2 ' 
M(r-2M) . 



Agg T 

:=A 4 

-r(r - 2M) , 

Ag&gth — A^ggdj = Ag^g = A^g^g = A, 



3r 2 



M 



sin z 9 . 



3(r- 2M) 
M 



3(r- 2M) 



sin 2 6 . 



r(3r - 2M) 



■ sin 



= -r(r - 2M) sin 4 9 - r 2 sin 2 9 cos 2 9 , 

A r g^ = A r ^gj, = A r ^g — Ag r ^ = A^ r g^ — A,p r ^g = Ag^^, 

= A^g rt j, = AtjftjyrQ = Ag^r = A^g^. = A^g r = r sin 9 COS 9 
R R -R M 

D tt\r — ~~ D tr\t — D rt\t — ~ i 
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B 



M 



rr\r 



(r-2M) 2 ' 



Bgg\ r = -Bg r \ e = B r g\g = T , 

— — -B^ r |0 = #r0|0 — r sin 2 9 , 
B(,M,h = -B^ = B^o = r 2 sin 9 cos ( 
M 2 (r - 2M) 

M(2r- 3M) 



B, 



tt\tt 



B tt\rr r 3( r _2M) ' 

_ M(r - 3M) 

^tr|tr " B tr \ rt - -^ {r _ 2M) 

M(r — M) 
B rt \tr ~ B rtlrt - r3(r _ 2M) , 



B 



M 2 



rr\tt 



r 3 (r-2M) 



B t 0\te — B t g\et — B et \ te — B 0t \ et — 



M(r - 2M) 
2r~ 2 



M(r-2M) 2 

OfMtA — DtMrht — t3 r l,t\trl, — Dm\m — r-^ Sill (7 , 



B 



t<t>\t<j> — n t<t>\<t>t — n 4>t\t<j> 
M(2r - M) 



2r 2 



rr\rr 



B r 



■\ee 



r(r - 2M) 3 
-1, 



B r 0\ r — B r g\g, r 



M 



Bor\rO — Bg 



\0r — 



sin 2 9 . 



r> i ) 



2{r-2M) 
2r-3M 
' 2(r - 2M) 

M 

~2{r-2M) 
2r-3M 
~2(r-2M) 



sin 2 9 , 



sin 2 9 , 



B 09 \ee = -r(r - 2M) , 

#00|00 = -r 2 cos 2 , 

Be4,\e4, = B e ^ e = -r(r - M) sin 2 6> , 

#00|00 = B<\>e\w = -r 2 cos 2 9 + Mr sin 2 9 , 

#00100 = -r 2 sin 2 6> , 

#00100 = ~r(r - 2M) sin 4 9 - r 2 sin 2 cos 2 

B r e\<f><fi = #0r|00 = —Bg,p\ r< fi 

= —B r ,h\j,g — B^Adrh = B^ r \j,g = —Ba 



Bg^r — B^g\ r ^ — B, 
= —Ba 



• \8r 



— — B r 
-r sin 9 cos f 



APPENDIX B: THE DIRECT PART OF ELECTROMAGNETIC AND GRAVITATIONAL SELF-FORCE 

In this appendix, we summarize the direct part of the vector and tensor field. The direct part of the field is obtained 
by integrating the direct part of the Green function 8 G^jj as same as the scalar case. 

a G$ } (x, ^) = ~ 0[V(x), x'] s u {A} (x, x')S(a(x, x')) , (Bl) 
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f VA(x,x>) (s = 0), 

i {A} (x,x') = I ^A(x,x')g ^{x,x') (s = 1), 

I y/ A(x, x')g^> (x, x')g vv i (x, x) [s = 2), 



(B2) 



where the suffix {^4} stands for the spacetime indices appropriate to the spin s of the field. 

From the above Green functions, we obtain the direct part of the field which is expanded around the particle 
location as 



= < 



cr- a (x, z(T TCt ))v a (T Ict ) 
g^ a (x,z(T ret ))v a (T iet ) 



0(y 2 ), 
-0(y 2 ), 



AGm 



(T-p{x,z{T Iet ))v^(T Ie t) 

g^ a (x, z(T rct ))gup(x, z(T ro t))w Q (T re t)w /3 (r re t) 



scalar 
vector 



(B3) 



0(y 2 ) , tensor , 



And then, we have 

F a [ s (l>fX } ](x) = cg a a (x,z cq )-^- \a- & {x,z cq ) + he 2 R a ^ s (z cq )v^ q a n (x, z cq )vf q } + 0(y) . 



e = J2a(x, z c 





q ueq 



1 



1 + 9-Ra / 375( z cq)Wc q CT ;/3 (a;, z cq )v2 q a.g(x, z cq ) + 0{y A ) , 



G 

where c and h depenend on the spin of the field as 

" (q 2 , 1/3) , 



scalar , 

(c , h) = { (-e 2 , ' - 2/3) , vector , 
(Gm 2 , - 11/3) , tensor . 



(B4) 
(B5) 

(B6) 
(B7) 



Here an exten sion of the four velocity v a (ro) necessary to define the projection tensor P a ^, as mentioned in the line 
following Eq. (1.3), is chosen such that 



V a (x)=g a a (x 1 z cq )v a (T cq (x)). 



(B8) 



It is n oted t ha t whe n wc consider the mode decomposition regularization for the self-force, the direct part calculated 
in Eqs. (1.26)^(4.29) is independent of spin. 



APPENDIX C: BASIC PROPERTIES OF THE MODE COEFFICIENTS 

In this appendix, we examine the general properties of the mod e coef ficients for the terms that appear in the local 
coordinate expansion of the direct force given in Eqs. ( 4.26| ) ~ ( 4.2S| ) . We show that the C and D-terms of the 
regularization couter terms vanish in the coincidence limit to the particle position. 

We first express £ in the form, 

e=eo+c 2 {<t>-<t> r )\ (ci) 

where £o is defined by 

£o = V2r (a-cos9) 1/2 . (C2) 
In terms of £o> R an d <fi, all the terms that contribute to the direct force F^ 11 will have the form, 



vy J J "° ^ Sr , (C3) 



R p (<t> - <j>')Hi r (R^Y 1 (f_\ n 

where we have replaced possible factors of the form 2fe in favor of polynomials in £q, R and <f>, and m, n, p, q and r 
are non-negative integers satisfying 
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to > , n > , l<p + g + 2r<3. 



(C4) 



This is because the highest order of the divergence in the direct force is quadratic, and we may focus only on terms 
of order up to O(y ) in the local coordinate expansion. 

Let us analyze the coincidence limit in detail. By using a modified version of Eq. (O) with £ 2 = but by taking 
account of 0{y 2 ) corrections, one can further reduce the above to the form, 



q+2ri£2q+2n+2r—2m—3jj2m+p 



(C5) 



Note that the 0(y 2 ) c orre ctions only change the original q to q + 2, hence it is enough to consider the above form. 
We can decompose ( C5 ) into a spherical harmonic series by using the formula, 



C, 



21+1 f 1 P t {n) 



dfj, = V2 (a - \/a 2 - l) 



(C6) 



This Cn is equal to the special case of C?(a) with p = —1/2 given by Eq. (D3). Introducing a variable z by 
e z — a + Va 2 — 1 (hence z > 0), Eq. ( |C(j ) is re-expressed as 



Ct{z) 



21+1 



_l -y/cosh z — ii 



d/i 



-Lz 



(C7) 



where L =£ + 1/2. Note that z oc \R\. 

As discussed in Section |y|, odd powers of the operator will give the harmonic coefficients which vanish after 
summing over m. Hence we only need to consider the case of q being an eve n no n-negative integer. Then it is 
straightforward to see that the operator in the harmonic decomposition of ( |C5| ) will give the factor L 2 in the 

coincidence limit of the angular coordinates followed by the sum over m. In general, 9?" will give rise to a polynomial 
of degree n in L 2 . Explicitly, we have 



2 T 



(C8) 



p=0 



where the factor Ap™' 1 is independent of L and An™" 1 = T(n + 1 / 2) / (^/tt T (n + 1)). The derivation of this formula is 
given in Appendix C of Ref. 



v(n) 



Thus we may replace (d^) q ^ 1 in (C5) by L 2j (1 < j < q/2 + n; note that q is even). With this replacement, let us 
consider the following 3 cases for the powers of £o separately: 

(1) 27V +1 :=2m + 3-(2g + 2r + 2n) > 1 (^ (2N+1) ; N> 1). 



In this case, to obtain the harmonic decomposion of £ 
to it: 



(2N+1) 



_! (C0shz-/i)( 2 »+ 1 )/ 2 



d/i oc 



, one simply applies [d/d(cosh z)] N 
d 



sinh zdz 



Taking account of the general form (C5), this will give rise to the harmonic coefficients as 



Lz 



\dj sinh zdz] N 
(C9) 

(CIO) 



where 1 < k < N. Since z ~ \R\ and 

2m + p = 2N + 2q + 2r + 2n - 2 + p 
>2N + q+(q + 2r+p)-2 
>2N + q-l>2N -1, 

the only term that remains in the limit R ~ ±z — > is the k — 1 term, and this implies 2m +p is odd. Since we 
know the leading divergence is quadratic, this will give a harmonic coefficient proportional to sign(i?) L which 
contributes to A,,. 
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(2) 2m + 3 - (2q + 2r + 2n) = 1 (Co" 1 )- 

In this case, the harmonic coefficents are non-vanishing only if 2m + p = 0. Since the result is independent of 
L, it contributes to B^, i.e., the linearly divergent term. 

(3) 2N + 1 := 2q + 2r + 2n - (2m + 3) > 1 ((% N+1 ; N > 1). 

In this case, since the harmonic coefficients of £,q N+1 will be finit e in the limit z — » 0, we must have m = p = 
which implies q + 2r = 2. Going back to the original form (|Cq), one can then redefine n by n + q/2 and the 
term of our interest takes the form, 



W 1-1 (»>!)■ 



(Cll) 



Using the formula (D4), we obtain the coefficients of the Legendre decomposition of £ ™ as 



t' 2 ' 
So 



i (L 2 -l 2 )(L 2 -2 2 )---(L 2 -n 2 ) ' 



(C12) 



where k„ = (-1)™ [(2n - 1)!!] 



and we have introduced the notation, 



to denote the L = £ + 1/2 mode coefficient of the Legendre expansion in the coincidence limit.^j Therefore, 
together with Eq. (C8), we have 



2n c 1n-\ 



°<h So 



(-1)"K„ 



(L 2 -l 2 )(i 2 -2 2 )...(L 2 -n 2 )^ o 



p ^ o (L 2 -l 2 )(L 2 -2 2 )...(L 2 -(n-p) 2 )' 



(C13) 



where f p is independent of L. We thus find the first term on the right hand side, which is inde pende nt of L, 
contributes to B^, while the rest seem to give non- vanishing contributions to D^. However, Eq. ( |C12| ) tells us 
that they are simply the Legendre coefficients of positive powers of £o, hence they vanish after the sum over £ 
is taken. 

In all of the above three cases, nothing contributes to C^. Thus, to summarize, we find the C and D-terms vanish 
and the A-term is proprotional to sign(i?). 



APPENDIX D: MATHEMATICAL FORMULAS FOR MODE DECOMPOSITION 

In this appendix, we give formulas neccesary for the harmonic decomposition of £ 2 ™~ 1 . For this purpose, we 
introduce a dimensionless version of £ by 

1= (a-cos£+^-0') 2 

Here, as defined in the text, 9 is the angle between {8, <fi} and {9', </>'}, and a > 1. In the following, we do not restrict 
n to be an integer. Our strategy is to consider a series expansion of £ by treating b as the expansion parameter. 
First, we consider the harmonic expansion of £q P where io = (a — cosf?) 1 / 2 , 

oo I 

U P = (o - cos §y = 2nJ2J2 G!(a)Y em (6, <^)F/ m (9', 0') . (Dl) 

£=Q m=-l 




When the m-sum is non-trivial, ■ ■ - \l is defined to be the coefficient after summating over m. 



1G 



The expansion coefficients C\ are obtained by the Legendre integration as 

C\{a) = J dfi(a-fiyP e ((j,), (D2) 

where Pe(fA) is the Legendre function of the first kind. As a > 1, we can expand (a — fj,) p in powers of p,. Then we 
obtain 

rv, 1 p-t r(fc - g/2 + l/2)r(fc - p/2 + £/2 + 1/2) 1 / 1 \ k 

WW- 2 p+ir(- P ) ^ r(fc + ^ + 3/2) A! U 2 / 

~ V r(-p)r(^ + 3/2) V 2 ' 2 '* + 2'a 2 

= a"*" 1 F -H h 1, - H h-,p + 2;l 

p+1 V 2a / V 2 2 2 2 2'^ ' a 2 

If 11*2*+ r(p+1)2 FT P I 1 + - -vl-- ) (I)'-!) 

+ l ° r(p-£+l)r(p + J?+2) ^ 2 + 2' 2 + 2 + 2' P ' 1 a 2 ' ' 



In the coincidence limit a — > 1 + 0, we have the two qualitatively different leading behaviors depending on the value 
of p as 

' -(a-l) p+1 , forp+l<0, 



P 1 r(p+1)2 (D4) 
(-lf2 p + l — LlP+i) forp+l>0. 

It should be noted that the divergent behavior persists in the mode coefficients for p + 1 < 0. 

Next, we consider the harmonic expansion of ((f) — 4>') n (a — cos8) p . We take into account only the leading order 
behavior in the coincidence limit. Hence we have the basic formula which converts a power of ((f) — <f>') to the 
(^-derivatives, 

(6 - (/>')( a - cos0) p = -^—^-{a - cos0) p+1 + 0(y 2p+3 ) . (D5) 
p + 1 d(f> 



Using this formula reccursively, we obtain 



[™/2] q„-2fc 



((f, - (f)T(a - cos§) p = £ 4 n) Sr^(« - co S 0)*+"- fc , (D6) 

where [n/2] denotes the maximum integer not exceeding n/2, and satisfies the recurrence relation, 

= l -T—- r4" ) -("- 2fc + 2 )4-i. (D7) 

k p + n-k + l ky ' " 1 v/ 

_(0)_/l, fc = 0; m „s 

-io, fe = l,2,3,.... (D8) 



This is solved to give 



(n) = r( P + i)r(n + i) (-i) fc 

lfc r(p + n-fc + l)r(n-2fc + l) 2 fc fc! 



(D9) 



Since each ^-derivative is interpreted as giving one factor of im to the harmonic coefficients in Eq. (Dl), we have the 
mode decomposition as 



6') n (a-cos0> = 27r£ £ C& p \a)Y tm (0, d>)Y t W A') , (D10) 



-i(n,p) t 

1=0 m —-l 
[n/2] 

= £ M«- 2 ^Cr- fc (a) > (DH) 

fc=0 
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With the formulas (D3) and (Dll) at hand, we can now write down the harmonic expansion of £ 2p to the leading 
order of the local expansion around (0', </)'). The result is 



E 



r(p + i) i (b 

T(p-n + l)n\ \2 



oo e 



27r E E 



0') (a-cos6») p - n 



(D12) 



where 



DL(-) - E 



r(p + i) i /6 



n=0 
oo n 

EE 

n=0 fc=0 



r a 



T(p-n + l) n! \,2, 

r(p + l)r(2n + 1) 



a) 



r(p + n - fc + l)r(2n - 2fc + 1) 2™+ fc n!fc! 



(-m 2 )"- fc Cj +n - fc (a) . 



(D13) 



The double sum in the last equation with respect to k and n may be simplified by introducing n :— n — k and summing 
over n and k, which can be now taken independently. Then we have 



oo / oo 



m m (a) = E E 



T(2n + 2fc + l) 1 f-b 



n=0 \fc=0 



r(p + i) 



r(n + /s + l) fc! V 4 / / r(p + n + l)r(2n+ 1) 



2, \ n 

—m^b\ 



C|+"(a) 



1 CXD 

= E 

1 _L 7, ' 



r(p+l) 1 / -m 2 b 



VT+b ^ T(p + n + 1) n! ^2(1 + b) 



(D14) 



Using Eq. (D4), the leading behavior of the mode coefficients (D14) in the coincidence limit a — > 1 + becomes 

1 1 



(a - 



T(p+l)T(p + n + l) 



1 / —m 2 b 



VTT5 ^r(p + n-l+l)r(p + n + f + 2)n!U + () 



if p+ 1 < 0; 
, if p + 1 > . 



(D15) 



Replacing p in the above by p — 1/2, we have the formula referred to in the text, Eq. (4.30). 
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